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ABSTRACT 


The present work Is a practical approach to study the 
nature of forces arising from one-dimensional sea-waves. The 
waves, both regular and random were generated by a wavemaker in a 
wave flume. The wavemaker is PC-controlled, Forces and moments 
acting on the vertical cylinder instal l^ed in the flume were 
measured by two properly mounted force transducers. The cylinder 
was vertically mounted as a cantilever with no support at the 
bottom. 

Both in-line and transverse forces acting on the 
cylinder have been measured. The transverse forces are generally 
less than in-line forces In case of bottom-supported structures. 
However, in the present study, the dynamic response of the 
cylinder and its vibrations due to the free end at the bottom is 
seen to result in higher values of transverse forces compared to 
in-line forces. Experimentally determined forces and moments have 
been compared with the forces and moments predicted by the Morison 
equation. 

Experimental results for regular waves are seen to match 
with results of linear theory. For random waves, water surface 
elevation, forces and moments have been recorded. Data has been 
analysed using Fast Fourier transform alogrithm to get wave, force 
and moment spectra. The theoretical spectra have also been 
computed using a linearisation of Morison equation. Experimental 
force and moment spectra are seen to follow the same trend as the 
theoretical spectra over a limited range of frequencies. 
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CHAPTER I 
i nlr oduct>l on 

The waves generated at the surface of seas induce forces 
on the structures installed. in them. It is of utmost interest for 
a structural designer to know the dynamic characteristics of the 
wave loading. 

The present work 'is a simple approach for studying the 
wave characteristics and the corresponding force and moment 
characteristics. The model is simplified in the sense that only 
one dimensional waves have been studied. The waves are generated 
in a wave flume by a wavemaker. Studies have been carried out 
with regular waves and random waves with Jonswap spectra. 
Theoretical and experimental force and moment spectra have been 
obtained and compared. 

In the present thesis, Chaper 2 deals with the 
literature review and the proposed theories available for regular 
as well as random waves. Chapter 3 deals with the experimental 
details and the procedure adopted for conducting experiments. 
Results have been presented in Chapter 4 and have been 
subsequently discussed. Chapter 5 presents conclusions of the 
present work. Appendix A contains computer programs written in 
Fortran 77 and the sea- language software for running the 
PC-controlled wavemaker. Appendix B deals with the calculations 
for regular waves. Some relevant definitions have been included 
in Appendix C. 
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CHAPTER II 

Literalure Review 


2. 1 Background Study 


Sea-waves can induce both in-line and transverse forces 
on an offshore structure. The lift forces are induced due to 
vortex shedding. The fundamental results of interest concern the 
inertia force, which is exerted on a body held in a uniformly 
accelerating ideal fluid and the drag force arising from an 
asymmetric distribution of pressure. The inertia force is written 
as , 


where. 


F . 
1 


F . 
1 


m 

P 

V 

du 

dt 


I, 

dt 


inertia force 
inertia coefficient 
density of fluid 
volume of body 
acceleration of fluid 


( 2 , 1 > 


A method of calculating wave forces on piles was first 
proposed by Morison, Johnson, O'Brien and Schaaf (1950) and this 
provided a landmark in wave force prediction methods. Morison and 
his co-workers suggested that the forces acting on a section of a 
pile due to wave motion is made of two components: a drag force, 
analogous to the drag on a body subjected to steady flow of a real 
fluid associated with wake formation behind the body; and an 
inertia force, analogous to tha.t on a body subjected to a 
uniformly accelerated flow of an ideal fluid. 
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Morlson equation Is thus expressed as follows: 

F = 0.5 pDC. UlUl + 0.25 prrD^ C ^ (2.2) 

where, 

D = diameter of cyl inder 

U = incident (undisturbed) velocity taken at the 

section centre 

= drag coefficient 

C = inertia coefficient 
m 

The application of Morison equation carries the implicit 
assumption that the body size is small relative to the wavelength 
so that the incident flow is virtually uniform in the vicinity of 
the body* When the body size is large, this equation is no longer 
val id. 

Hydrodynamic coefficients for forces on a circular 

cylindrical structure have received considerable attention, A 

systematic approach in this direction was made by Keulegan and 

Carpenter (1958) who tested horizontal tubes at the node of a. 

standing wave in a tank so that the tubes essentially experienced 

one dimensional flow. It was observed that C, and C depend on a 

d m 

period parameter known as Keul egan-Carpenter number defined as 
U T 

^ , where U is the maximum water particle velocity, T = the 
D m 

time period and D = diameter of the cylinder 

Relatively recent is the finding of Sarpkaya (1976). In 

his study, tests were conducted on horizontal cylinder in a 

vertical U-tube in which water was harmonically oscillated. He 

also found C_, and C to be dependent on the time period, 
dm 

Sarpkaya (1974) and Tutor (1974) have measured the lift 
(transverse) forces on their test cylinders. They found the 



presence of several lift frequencies in a. given wave cycle. A 

dominant frequency of vortex-shedding was shown to exist at twice 

U T 

the wave frequency. When Cj^ was plotted against — the maximum 
lift force was observed to be as high as 1.3 times the maximum 
in- 1 ine drag force. 

Flexible oscillating cylinders have been studied by 
Laird (1962J. It was found that dynamic response of the body 
could lead to high lift, as high as 4,5 times the drag on a fixed 
cylinder in an uniform flow. Thus the degree of rigidity and 
support conditions of the structure control significantly the lift 
character is t ics . 


Isaacson (1974) observed that for ■ |^ ■ < 2, the flow did 

U T 

not separate. For ^ ^ 4, a pair of eddies formed, which 

attached to the cylinder and moved symmetrically back and forth 

U T 

during each half cycle of wave. For ' ■ in the range 6-15, one 
eddy separated completely during each half cycle. Isaacson 

calculated the first five components of the lift coefficients. 
These were written as a function of the KeuI egan-Car penter number 



These findings gave impetus to research in wave 
mechanics and form the main background for the present work. 


a. 2 Aim of Study 

The aim of the present study is to measure the wave 
forces for the case of regular and random waves and to examine 
whether these forces can be predicted, 

2. 2. 1 Regular waves 

When a cylinder is subjected to one dimensional regular 



wavas of various frequancies and small wavehaights, it aKparlences 
inline and transverse (lift) forces. Maximum In-line force and 
moment can be theoretically predicted using the following 
important variables: 

= dimensionless wave steepness 

= dimensionless still water depth 

= relative pile roughness 
= a form of Reynolds Number. 

A appropriate wave theory is assumed to determine the 

theoretical value of the hydrodynamic coefficients C. and C . 

d m 

From the calculations presented by Dean < 1965a), curves have been 

C D 

prepared in the form of isolines for various ratios. These 

plots are available in the form of wave steepness as a function of 
the dimensionless depth ^ ■ ^ . Thus one can get theoretical maximum 

tr 

in-line forces and moments acting on the circular cylinder 
subjected to a unidirectional regular wave. 

Sarpkaya ( t976b, 1977a) has measured the transverse 

forces acting on smooth and sand-roughened cylinders for a wide 
range of Reynold numbers, Keul egan-Carpenter numbers and relative 
roughnesses. These results are available in graphical form. They 
have been used for comparison in the present work. 

2. 2. 2 Wave forces due to Random Waves 

Ocean waves are random in nature. Thus getting 
statistical characteristics of random waves is of concern in ocean 


H 



d 

gT^ 


k' 

D 

HD 

Tt> 


engineering. 


Random waves result in random loads on offshore 



structures. A random wave can be decomposed into several regular 
waves each of unique amplitude and frequency. 

In the present thesis, random (JonswapI wave spectra 
have been dealt with and the consequent loading namely, force and 
moment spectra have been determined. 

Before going into details of spectral analysis, it Is 
worthwhile to discuss some of the statistical properties of random 
waves. These are given below 

Variance, a ^ - IS (f) df 

K J X 

where S <f) is a continuous function of frequency. S (f) 

X X 

represents the contribution to variance (or "energy") of signal 
x<t) due to its content within the frequency range f to f + df. 
Symbolically we can write as 

f +df 

- 2 I l*nl^ '2-3' 

n 

where random wave signal can be expressed in a Fouries series as 
fol lows 

x<t) = Re ^ (2.4) 

After the signal is recorded we can easily get the spectra using 
Fast Fourier transform (FFT) alogrithm. 

Similar spectra can also be derived for the force and 
moment signals recorded in the experiments. 

One more important characteristic of random waves is H 

S f 

the significant wave height. It is a value above which ^ of the 
largest wave components are present. 

In order to find the nature of the probability 
distribution for the wave heights, we define a spectral width 


parmeter 



7 


€ * ll - 


m m, 
o 4 


(2.5) 


where 


■o ' / 


= J f^.S <f) df 


= J f^.S (f) df 


A value of € = 1.0 indicates Gaussian Distribution and 
indicates Rayleigh Distribution. 


2. 2. 3 Alternative approach lo represent one dim&nslorial spactra. 

The free surface elevation can be written in general as 


foil ows 


Z A„ .KP {i[k„ 


Cosie »x + K Sin (6 )y - 2nf 
n n n 




( 2 , 6 ) 


where = complex amplitude. 


Generalised three dimensional spectrum can be written as 

S(K,f,©); such that S (K, f ,S ) dK. df . d© represents the contribution 

2 

to the variance due to component waves with wave numbers 

between K and K+dK, frequencies between f and f + df and 
directions between & and & + dd. 

Thus we have 


K +dK f +df e *de 
n n n 


S(K, f ,d ) dK df dd 


5 5 


<X) to n 

VAriance o* ^ = J J J S ( K, f ,(9 ) dK df d© 

o o -fj 


(2. 7 ) 


This general three dimensional spectrum is too unwieldy 


for practical use. So, in two dimensions it can be written as 



foil ows 


2 00 n 

” J J S{f,d) df ae (2.8) 

O “fl 

where the two dimensional directional frequency spectrum S(f,d) 
can be obtained from the following numerical integration 

06 

S(f,^) = J S{K,f,^) dK (2,9) 

o 

The directional frequency spectrum S<f,0) has found 
general application in engineering. Similarly, other spectra such 
as directional wave number spectrum S<K,6) can also be obtained. 


Now from these 

two dimensional 

spectra. 

one 

d imens ional 

wavenumber and frequency 

spectra may be 

obtained 

by 

integrating 


the corresponding directional spectra over d, giving 


n 

S(f) = J S(f,e) 66 (2.10) 

-rr 

n 

S(K) = J SCK,©) 66 (2. 11) 

-rr 

2 

and the variance is given as 

n 


frequency 


O' = f S(f) df 
n J 

o 

The relationship 
spectra is given 


oo 

= ; 

O 

between 

as 


S(K). dK 

the directional 


wavenumber 


and 


s(K,e) = sif,e) (2.12) 

Zn G 

where 

C- = group velocity 
G 



Transformation of random wav© train due to refraction 
was described originally by Longuet-H igg ins (1956, 1957b). It has 



been shown that the energy flux between neighbouring orthogonals 


of any component wave train remains constant. Hence, 

S{f,d) df de C_b = S (f ,0 ) df de C b 
G QOO ooGo 

o 

This is equivalent to assuming linearity and hence 
transfer between any pair of frequencies. It can also 
that 


<2. 13) 

no energy 
be shown 


bK dS = b K d© = constant (2.1A) 

o o o 

Substituting Equation (2.1A) into Equation (2.13) and bearing in 
K ^o 

mind that p— = -pr- (since f is constant), we obtain 

K G 

O 

c c s(f,d) = c s (f, © ) 

G o G o o 

o 

= Constant along an orthogonal. 

In the present work random waves have been generated 
using Jonswap spectra with a randomized initial phase. The 
Jonswap spectrum is basically a modification of the 
P lerson-Moskowl tz spectrum (196A). The P ierson-Moskowitz spectrum 
can be written as 


SI f ) 



(2. 15) 


where a and ft are constants and f is frequency. The constant ft 
depends on velocity U, and is valid for a fully developed 
condition. 


The Jonswap spectrum accounts for fetch- restr ict ions and 
provides a much sharper peaked spectrum. Jonswap spectrum is 
given as 


ag 


{2fT) f ® 



S<f ) 


exp 


(2. 16) 



IV 


where a = exp 



Here, f is the 
o 

ratio of maximum 
waximum spectral 
Based 

the theoretical 


I - ( 

= 0.07 for f < f 
a o 

. = 0,09 for f > f 

b o 

peak frequency at which S{f) is maximum and 
spectral density in the Jonswap spectrum 
density in the Pierson-Moskow i tz spectrum, 
on these spectra, we are in a position to 
spectra of random wave forces. 


y is 
to the 

discuss 


2. 2. 4 Random wave for ces 

So far we have discussed the statistical nature of waves 
themselves. We turn now to consider the consequent loading on a 
structural element. The basic assumption here is that the Morison 
Equation (2,2) is always valid. Therefore a statistical 
description of fluid velocities and accelerations is required. 

Using linear theory, velocity and acceleration of wave 
trains are given as 


U(t) = H (f) T)(t) 
u 

a<t) = Hg{f) T 7 <t) 
where rjlt) is the input wave height signal. Further, 


(2, 17) 


H (f) 
u 


Cos h IKlZ+d) 
Sin h (Kd) 


H <f ) 
a 


- i4n^f 


Cos h {K(Z+d) 
Sin h (Kd) 


Consequently the corresponding velocity and acceleration spectra 


are given as 


(2, 18) 
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F i g . 2 . 1 


Flow Chart to Calculate Velocity Spectrum. 
Acceleration Spectrum and Force Spectrum 








Now from the above probabilistic properties of velocity and 


acceleration spectra, Morison equation gives the statistical model 
of force as follows 

Fit) = uit) (u(t)| 4 ait) {2.19) 

Coefficients and K are drag and inertia terms 

respectively and are assumed as constant for the whole depth of 
the cylinder. Here, 



and p, D, C. and C are density of fluid, diameter of cylinder, 
d m 

drag coefficient and inertia coefficient, respectively. 

The spectral density of forces can now be found using 

the flow chart given in Figure 2,1. 

In the figure given above, Sp i f ) is the spectral density 

2 

of predicted force on a section of the cylinder. Parameter <7^ is 
the area of the velocity spectrum at a fixed elevation 2. 

In the present thesis total force acting on the cylinder 
has been measured. The experimental force spectra hence 

correspond to the total force spectra. 

In order to get the theoretical spectrum for the full 
cylinder the above predicted sectional force spectra given in 
Figure 2.1 should be integrated numerically over the entire depth 
of immersion. 

Borgman il967,b) has provided force as well as the 
moment spectra for a single cylinder extending from the sea-bed 
and piercing the free surface. They are as follows 



SFT<f) 


S ( f ) 
T7 


{ _ j- 2nfK o 

n [ SU TKi - m J/u‘2) Cosh <K<Z.d))d2 J 


4 .. 2 


< 2 . 20 ) 




S { f ) 
ri 


(2nf) K. ■] 

{ _ IT 2nfK . o 1 

- rnT^TT r <Z + d)<y <Z) Cosh<K <Z + d ) ) dZ 

rr[^Sinh(Kd)J^ u J 


(Zrtf) K^[(Kd) SinhCKd) + 1 - Cosh(Kd)] 
Sinh (Kd) 


< 2 . 21 ) 


Here S <f) is the known experimental wave spectrum. 
7J 

S^„{f) and Su_{f) are the total force and moment spectra. cr <z) 
can be derived from the predicted velocity spectra at each section 
of the cylinder. Other constants and variables have been defined 
ear I ier , 

The goal of the present work is to derive the 
theoretical spectra for total force and moment using Equations 
(2.20) and (2.21) and to compare them with the experimentally 
measured profiles. 



CHAPTER III 


EXPERIMENTAL DETAILS 

3* 1 Equi pwient 

Experiments were conducted in a wave flume having 24m 
length, 45 cm width and 90 cm depth and having glass walls on the 
two sides* All experiments were conducted at a still water depth 
of 70 cm* The wavemaker consisting of control box, interface box 
and wavegauge was installed at the upstream end of the flume so 
that waves could be generated at the desired parameters set 
appropriately* An aluminium cylinder of 5 cm diameter was 
fabricated and installed in the flume as shown in the Figure 3*1* 
It was fitted with two transducers vertically separated by a 
distance of 5 cm* The idea behind putting two transducers is to 
measure forces as well as moments simultaneously* A well 
supported top end also facilitates supporting the cylinder in the 
experiment* The wavemaker and the cylinder in mounted position 
can be seen in photographs Nos. 1 to 5* 

The measurements were taken for a given wave signal 
input forces in longitudinal (in-ltne force) and transverse (lift 
force) directions* 

3. 2 Instrumentation 

The experimental setup to determine the wave 
characteristics is shown in Figure 3*1 and Figure 3*2. The 
Instrumentation consists of wavegauge, control box, wavemaker, two 
force transducers, interface box and oscilloscope. Their 
specifications have been presented in Section 3.4, 
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FIG-3-2 LONGITUDINAL SECTION OF THE CHANNEL AND 
PLAN VIEW OF THE CYLINDER INSTALLED 
(Dimensions ore not to the scale) 
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Photo, 1 



Photograph showing the wave flume 
with piping for filling and emptying 
the wave flume and wave monitor 
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and control box 



Photo. ? Photograph showing the storage 
oscilloscope and PC 



Photo . 


Photograph showing the wavemaker 


Photo. 5 


Photograph showing the cylinder mounted 
in position with the force transducers, 
channel selector, conditioning amplifier 
and wave probe 
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A wavegauge is used to measure the water surface 
elevation (t)) at any instant (t). It is of capacitance type and 
is linear in the working range with a ± 10 volt output. In order 
to convert the data obtained (in volts) to cm units, the wavegauge 
is calibrated using an oscilloscope. 

The control box contains the power amplifier for the 
wavemaker, along with signal conditioning, force feedback, wave 
absorption and nonlinearity correction circuitry. 

The wavemaker consists of a "flap* which is a wedge 
shaped paddle pivoted about its tapered end on a fabric hinge 
fixed to the bottom leading edge of the wavemaker box. The paddle 
is mounted in the box and extends the full width of the tank. The 
front of the box is sealed with a waterproof fabric gusset, behind 
which the paddle moves. The paddle Is driven by an electric 
servo-motor via a ber r i I 1 iura-copper drive belt which runs over a 
sector attached to the top of paddle. 

The two force transducers are of model 8200 (Bruel and 
Kjaer, Denmark) and sense forces in the range (-1000 N to 5000 N). 
Their specifications have been given In section 3.4. 

The interface box acts as the junction between the 
wavemaker and the outside world. All the signals to and fro the 
interface card in the back of PC are run by ribbon cable to the 
interface box. The box also contains a sine wave generator for 
running the wavemaker without using the PC. The system is run 
from an internal power supply. 

The interface card is manufactured Advantech (Model PCL 
812) and plugs into the expansion slot of any IBM-compatible PC. 
It provides two 12 bit analogue output channels, 16 analogue 
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inputs, 16 digital inputs and 16 digital outputs. 

The instrumentation can be seen in photographs Nos. 

1 to 5* 

3. 3 Procedxire 

The 5 cm diameter aluminium cylinder of height 90 cm was 
installed in the flume as shown in the Fig. 3.2. The cylinder was 
connected to the force transducers (8200 type) in the in-line as 
well as transverse directions in separate experiments to measure 
in-line forces and lift forces, respectively. 

The wavemaker operates for 1.5 minutes for each run. 
Data are recorded at 0.0625 seconds intervals. 

For regular waves of required amplitude and frequency, 
the wavemaker was run with suitable software through a PC. The 
data output (wave gauge signal and force signal) were recorded 
simultaneously. A given input signal was run twice to get the 
in-line and transverse (lift) forces. 

Similarly, for random waves the wavemaker was run to 
produce a Jonswap spectrum. All time-based signals were analysed 
using the Fast Fourier transform alogrithm; a computer program for 
this alogrithm is given in the Appendix A. 

For force signal analysis, the transducers were 
connected to the conditioning amplifier for amplification of 
signals and conditioning amplifier was connected directly to the 
Interface card to store the data. 

A total of eleven experiments have been conducted in the 
flume. Out of these, six are for regular waves at different 
amplitudes and frequencies, and five experiments have been 
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conducted for Jonswap spectra 
frequencies. The parameters for 
tabulated in Tables 1 and 2 
respective! y . 

The results have been 
graphical form, 

3* 4 Specifications 
3* 4* 1 Qsci 1 loscope 
Model number 

Manuf ac turer 
No, of input channels 
I nput coupl lag 
Vertical sensitivity 

Vertical gain accuracy (dc) 
Time base accuracy 
Time base range 
Maximum sampling rate 

3, 4* 2 Conditioning amplifier- 
Model 

Manufacturer 
Maximum sensitivity 
Output modes 

Frequency limits 


, set at five different peak 
the different experiments are 
for regular and random waves, 

produced in tabular as wel I as in 


: HP 54501 A 

Diziting Oscill oscope 

: Hewlwett Packard 

: 4 

j AC, DC 

: Maximum 5V/div 

Minimum SmV/div 

: +1.5% 

: 0.005% 

: 2 ns/dv to 5 s/div 

10 MSa/s 

: Type 2626 

; Bruel and Kjaer, Denmark 

: 1000 mV/ PC 

s Direct and transformer 
coupled. 

: Switchable low and high 

frequency. 
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Calibrated output rating 
Maximum output 
Level indicator 

Rise time 
Recovery time 
Environmental condition 
Temperature range 
Humidity 

Power requirements 

3. 4. 3 Force Transducer 
Type 

Manufacturer 
Force range 

Linear i ty 

Charge sensitivity 

S 1 if f ness 

Weight 

3. 4. 4 PC 

Type 

Manufacturer 


0.001-10 V/unit selectable 
in 20 dB steps. 

10 V ( 10 mA ) peak from 
0.3 Hz to 30 KHz. 

"overload* Led lights when 
input or output levels 
exceed lOV peak. 

~ 2,5 V/^S 

< 200 (jS 

-10 to + 55"C 

0 to 90% RH 

100- 240V (50-400 Hz t + 10% 
AC.7VA complies with safety 
class 1 of EC 348. 

B & K Type 8200 
Bruel it KJaer, Denmark 

1 KN tensile to 5 KN compre- 
ssive. 

< ± 1% of maximum force 

4.0 PC/N 

5.0 X 10® N/m 
21 gram 

386 SX 

Wipro (Pvt) Ltd., Bangalore 
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3. 4. 5 Interface Card 
Model 

D i g i t ieat ion 
No. of channels 
Maximum sampling rate 


: PCL 812 

; 12 bit 

: 16 A/D 

: 5 KHz 


3- 4. G Wavegauge 
Type 

Output voltage 
Manufacturer 


; Capacitance 
: ± 10 vol ts . 

t Armf ie I d , U . K, 
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CHAPTER IV 

RESULTS AND DISCUSSIONS 

The results have been presented in tabular form In 
Tables 1-9, and as graphical plots in Figures 4. 1 - 4.8 at the end 

of the chapter. 

Regular waves - 

In order to calculate the drag coefficient (Cd), the 

inertia, coefficient (C ) and the lift coefficient (C, ), Reynolds 

tn L 

number and Keul egan-Carpenter number for each experiment are 

required. These are listed in Table 3. Table 4 shows the 

comparison between maximum in-line force and maximum in-line 

moment measured and those predicted by theory. Sample calculation 

for determining the maximum in-line force and moment are given in 

Appendix B. Table 5 shows a comparison between the lift force 

found experimentally and that predicted by linear theory. It is 

observed that the lift force as well as in-line force increase 

with increasing frequency and waveheight. 

In the above analysis the roughness parameter k' has 

been assumed to be very small since the cylinder material is 

aluminium. From the shore protection mannual (vol II), C , , C and 

a in 

C, are found to be 1.2, 2.0 and 3.3 respectively for all the six 
L 

experiments for regular waves. 

1 

The recorded signals for the regular waves have been 
presented in the Figure 4.1. Visually it can be observed that 
force signals are not in phase with the water surface elevation 


signal 
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A comparison of Tabls 4 and Table 5 shows that the lift 
force is more than the in-line force. This is due to the dynamic 
response of the cylinder which is not bottom-supported. 

Random waves ; 

In Tables 6 and 7, rms values of forces and moments have 
been presented. It is observed that out of five peak frequencies 
chosen, the maximum rms value of in-line force and moment is 
attained at a peak frequency <fo) of 0.85 Hz. The lift force and 
moment decrease continuously with a drop in the magnitude of the 
peak frequency. 

The spectral width parameter (el has been computed and 
tabulated In Tables 8 and 9. These are seen to be closer to unity 
and hence represent a near-Gauss ian probability distribution 
function. 

The recorded signals for random waves have been 
presented in Figure 4.2. Figure 4.3 shows the wave spectra for 
water surface elevation. Predicted and experimental spectra of 
the in-line forces and moments have been compared. This is shown 
in Figure 4.4 and Figure 4.5, respectively. The predicted spectra 
for forces and moments have been obtained from the formulas 
presented in section 2.2.4. Figures 4.6 and 4.7 show the spectra 
for forces and moments in the transverse (lift! direction, 
respectively. Lastly, Figure 4,8 represents the five harmonic 
components for drag and lift force and their phase difference with 
respect to water surface elevation. 

A comparison between experimental in-line force spectrum 
and that predicted by linear theory can be seen in Figure 4.4. It 
can be Inferred that at higher peak frequency, theoretical spectra 





are close to the experimental spectra. In this present approach 
for the calculation of theoretical spectra, the values of drag 
constant and inertia constant (K^ ) have been assumed as 
constant at all depths of water in the flume. Despite this 
assumption there is a good agreement between theory and 
exper iments. 



Similar calculations 

have been carried 

out 

for 

the 

in- 1 ine 

moment in Figure 4.5. 

The same trend as in 

the 

case 

of 

In-line 

force is observed here. 

Spectral analysis 

shows 

that 

all 


the peaks of theoretical and experimental force and moment spectra, 
are around the same respective frequencies. 

Predictions for the lift force and moment spectra in 
case of random wave are not available. The experimentally 
obtained spectra for lift force and corresponding moment have been 
presented in Figures 4.6 and 4.7, respectively. As before, the 
peak of the force and moment spectra is around the same frequency 
as it is for corresponding water surface elevation spectra. 

Figure 4.6 shows that the input water surface elevation 
signals and the corresponding force signals are not in phase and 
their harmonic components have phase differences. 
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Table 1 


Parameters for Regular Wave Experiments 


Exper iment 

Ko, 

F requency 

Hz 

Uavehe 1 ght 

Cm 

1 

l.OO 

8.00 

2 

1,00 

12,00 

3 

t.25 

8,00 

4 

1. 25 

12.00 

5 

0. 75 

8.00 


6 


0. 75 


12,00 



30 


Table g 

Parameters for Random wave CJonswap spectral Experiments 


S(f ) = 


_ exD f- 5 f ' 

1 

( 2*7 ) 

5 1 « 1 *o •» 

J ^ 

a *= 

exp 

1 

0 

ro 

0 




0.07 for f < f 


ct = 

o 



l-b • 

0.09 for f > f 

o 


f 

o 

peak frequency (Hz) 


a - 

0.0081 



g 

9*81 m/s^ 


Experiment 

Peak frequency 

r 

No. 


f (Hz) 
o 


7 


0.75 

5.00 

8 


0.85 

5. 00 

9 


1.00 

5.00 

10 


1. 10 

5.00 

11 


1.25 

5.00 
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Table 3 


Reynolds numbers and lCetjQ.e9an-Carpent.er nunbers 


Expt. No. 


maximum 

surface 

velocity 

U 


Reynolds No. Keul egan-Carpentei 

No. 

U D U T 

Re = K* = — 


cm/s 


1.00 

9. 47 

29.74 

14872 

5.95 

1.00 

12, 60 

39. 60 

19803 

7,92 

1,25 

8.33 

32.71 

16356 

5.23 

1,25 

10,89 

42. 76 

21380 

6,84 

0.75 

8. 42 

19.83 

9919 

5.29 


6 


0.75 13,23 


31. 17 


15585 


6,23 
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Table 4 

Experl menial and predlcled in— line force arid momenl for 
regular wave experiments 


Expt. f 

No. 

(Hz) 

Measured 

Wave height 
at SWL 

Maximum In- 

line force 

Maximum moment in 
In-line direction 


cm 

Exper imen- 
tal force 
(N) 

Predicted 

force 

iH) 

Exper imen- 
tal moment 
(Nm) 

Predicted 

moment 

(Nml 

1 1.00 

9.47 

2.38 

1.86 

1.68 

0. 96* 

2 1.00 

12.60 

3.85 

2.56 

2.88 

1. 39 

3 1. 25 

8.33 

3.55 

1.61 

2. 99 

0. 94 

4 1. 25 

10.89 

4.85 

2. 17 

4.23 

1. 30 

5 0.75 

8.42 

2.08 

1.57 

1. 39 

0. 70 

6 0. 75 

13.23 

4. 21 

2.54 

3. 35 

1. 20 


6 
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Table 5 


Experimental and predicted lift force and mciment for 
regiilar waveexperiments 


Expt, 

No. 

f 

(Hz) 

Wave height Maximum Lift force 

at SWL 
< cm) 

Maximum moment i 
Lift direction 




Exper imen- 
tal force 
(N) 

Predicted 

force 

(N) 

Experimental 

moment 

(Nm) 

1 

1.00 

8,42 

3.59 

5.20 

3.50 

2 

1. 00 

13.09 

6.90 

9.05 

6.74 

3 

1. 25 

8.46 

3. 94 

6. 18 

3.82 

4 

1, 25 

11.21 

7.64 

10,56 

7.34 

5 

0.75 

8. 74 

3.06 

2.27 

3.00 

7 

0.75 

12. 89 

6.88 

5.61 

6.78 
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Table 6 

Experimental rms value of In-line force and moment for 

random wave experiments 


Expt . 
No, 

Peak frequency 
f 

o 

(Hz) 

Rms value of 
water surface 
elevation 
measured at 

SWL (cm) 

Rras value of 
In-line 
f orce 
(N ) 

Rms value of 

I n- 1 ine moment 
(Nm) 

7 

0. 75 

2. 34 

0. 42 

0.56 

8 

0. 85 

1, 84 

0, 98 

0, 83 

9 

1.00 

1. 33 

0.50 

0.32 

10 

1. 10 

1.11 

0.41 

0, 26 

11 

1.25 

0,84 

0.31 

0.20 
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Table 7 


Experimental rms value of Lift force and moment for 
random waves experiments 


Expt . 
No- 


Peak frequency Rms value of Rms value of 
f water surface lift force 

® Cl X • 

E I e V at ion 

(Hz) measured at (N) 

SWL (cm) 


Rms value of 
lift moment 

(Nm) 


7 

0.75 

2. 31 

1.55 

1.51 

8 

0.85 

1.81 

1, 26 

1.23 

9 

1.00 

1. 32 

0.86 

0.84 

10 

1. 10 

1,11 

0,61 

0,59 

11 

1.25 

0.83 

0.51 

0. 49 
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Table 8 



Spectral width parameters for 

In-line force 

and laonient 

Expt. 

No. 

Water surface Elevation 

Force 

Moment 

€ 

e 

e 

7 

0. 77 

0. 79 

0. 78 

6 

0, 72 

0.65 

0.65 

9 

0. 65 

0. 85 

0.86 

10 

0, 60 

0. 84 

0.84 

tl 

0.57 

0.87 

0.87 
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Table 9 


Spectral width parameters for Lift force and moment 


EKpt. 

No, 

Water surface Elevation 

Lift Force 

Moment in transverse 
direction 


€ 

€ 

€ 

7 

0.73 

0,92 

0, 92 

6 

0.70 

0. 90 

0.90 

9 

0,62 

0. 90 

0. 90 

10 

0,60 

0, 92 

0. 92 


11 


0.56 


0.89 


0. 89 





Measured Signals of Water Surface Elevation 
Coefficient and In-line Force for Expt . 1 

(Regular Waves) 




-015 



Measured Signals of Water Surface Elevation 
Coefficient and In-line Force for Expt. 2 
(Regular Waves) 
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Measured Signals of Water Surface tieva 
Coefficient and In-line Force for Expt . 
(Regular Waves) 

















Measured Signals of Water Surface Elevation 
Coefficient and In-line Force for Expt . 6 
(Regular Waves) 






EiH 




asured Ra 










Measured Random Wave Signals of Water Surface 
Elevation. Lift Force and In-line Force for Expt . 9 
(f„ = 1.00 Hz) 
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Measured Random Wave Signals of Water Surface 
Elevation, Lift Force and In-line Force for Expt . 11 

<f„ =1.25 Hz) 
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A Comparison between Experimental and Predicted In-line 
Force Spectra for Expt . 7 (f^ = 0.75 Hz) 

(Dotted line represents the predicted curve) 
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Fig. 4.4.d A Comparison between Experimental and Predicted In-line 
Force Spectra for Expt . 10 (f^ = 1.10 Hz) 

(Dotted line represents the predicted curve) 
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Fig. 4 . 8 . a Five Harmonic Components of Force Spectra 
f « 0.75 Hz tExpt. 7> 
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Fig. 4.8.b Five Harmonic Components of Force Spectra 
f = 0.85 Hr <Expt . 8) 
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Fig. 4 .8.C Five Harmonic Components of Force Spectra 
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4-8. d Five Harmonic Components of Force Spectra 
f » 1 . t 0 Hz <Expt .10) 

NOTE: NO. WITHIN BRACKETS REPRESENTS PHASE ANGLE 
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Fig. ^.8.e Five Harmonic Components of Force Spectra 
f « 1 .£5 Hz (Expt. 11 ) 
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CHAPTER V 

CONCLUSIONS 

Based on the results obtained, the following conclusions 

have been arrived at, 

(1) The maximum in-line force and maximum in-line moment 

increase with increasing frequency of regular waves. 

(2) The maximum in-line force and maximum in-line moment 

increase with increasing waveheight of regular waves. 

(3) The lift (transverse) force may be larger than the in-line 
force if the cylinder is free to move at one end. Hence 
in design, lift forces can not be ignored even in one 
dimensional flow fields. 

(4) The experimental in-line force and moment spectra match 
with the predictions of Borgman( 1967, b) at higher peak 
f requenc i es . 

(5) With increasing peak frequency, water surface elevation 

spectrum becomes increasingly sharp, whereas the 

corresponding force and moment spectra remain nearly 
Gaussian and maintain an almost constant value of spectral 
width parameter (iSi 0.9). 

(6) In case of the random waves, the harmonic components of 
forces are not in phase with their corresponding water 
surface elevation components. 
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APPENDIX - A 


LISTING OF PROGRAMMES 


PROGRAM FFTAC 

CALCULATION OF POWER SPECTRUM OF SIGNALS USING FFT. 

dimension a<?200> .power <5200) .auto <5200) .alag<5200) 
dimension b<5200).a2 <5200). time <5200) 
complex x<5200) 
common / tr /troom 


print *, 'what are nu ndata del* 
read «. nu. ndata. del 
scale=10 . 0 

n=2**nu 
tt=del#<n-l) 
pi=5. 1415926 

open<un it=21.file='p'.status='old’) 
read<21 .*) <time<i).a<i},isl .ndata) 
close<uni t=21 .status keep ' > 

remove stray mean. 

call trap<a. ndata. e) 
do 65 i =1 . ndata 
a< i )=a< i )-e 
cont i nue 

filter data <Tukey- Hanning formula). 

ttl=0.5*del*<ndata-l ) 
do 75 i=l. ndata 
tl=del*<i-l) 
wt = l . 0 

a2< i )=a< i )**2 
a( i )=a< i )*wt«scale 
cont i nue 

call trap<a2 .ndata .variance) 

do 45 i=ndata+l.n 
a< i )=0 . 0 
cont i nue 
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call fft<a.b,x,nu,del > 

do 100 i *1 ,n 
r«real <x< i ) ) 
sc innag(x( i > ) 

power ( i >*(Bqrt(r#*2+s*#2> >/tt 
100 continue 

print *, ’power spectrum calculated* 
delf=l .0/tt 

open < uni t =2 2 , f i le= ' power . dat ’ , status = * unknown ' > 
fnyquist=l .0/(2 .0«del ) 
do 8> i^l .ndata 
i f <b( i > . ge . fnyqui st > go to 86 

85 continue 

86 i 

write (22.1130) 

1130 format(/2x, ’ frequency in Hz ; amplitude dimensionless'.//) 
write(22.1100) 

1100 format (3x. ' frequency* , 1 lx . ' ampl i tude * , //) 
write (22. 11 50) (b( i ) , power( i > . i=l . i j) 

1150 format(3x.el2.5.8x.el2.5) 

write(22.1190) fnyquist 
1190 format (/ 5x. 'Nyquist f requency= ’ . el 2 . 5 ) 

close(unit=22,status='keep’ ) 

stop 

end 

subroutine f ft(a.b.x. nu.del ) 

dimension x(3200>,b(3200).a(3200) 

complex x.u.w.t 

n=2^^*nu 

pi=5. 141592653 

tt=(n-l)*del 

do 5 j j=l .n 

b( j j) = ( j j-l)/tt 

5 continue 

do 6 i=l.n 

x( i ) =cmplx(a( i ) . 0 . 0 ) 

6 continue 

do 20 l=l.nu 
le=2**(nu+l-l ) 
lel=le/2 
u=(1..0.) 

w=cmpl x(cos(pi/fl oat (lei)) ,-sin(pi/float( lei ))) 


do 20 j=l.lel 
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c PROGRAMME TO CALCULATE WAVENUMBER & CELERITY 

C RANDOM WAVES 

c CALCULATION OF PHASE DIFFERENCE BETWEEN TWO SIGNALS. 

dimension si C 1 100) .s2< 1 J 00 ) , b( 1 1 00) . powerl ( 11 00> , power2 < 1 1 00) 
dimension al 1 ( 1 1 00) ,s22 < 1 1 00) 
complex xl ( 1 1 00) , x2< 1 1 00) 

dimension pl2 ( 1 1 00) . t ime( 1 100) .pcrdlOO) .pH ( 1 100) .p2 2(1 100) 
dimension akdlOO) .cel (1100) .phased (11 00) 

pi*=J.1415 

nu=10 

n=2*»nu 

pi=?.14H926 

open(uni t=6 . f i le= 'ph .dat ' . status^ *old' > 
read(8.«) ( t ime( i ) ,al 1 ( i ) . i =1 .n) 
read(8.*) (time( i ) .a22( i ) . i=l .n) 
seal e=l . 0 
do 20 i=l.n 
al(i)=all(i)«scale 
a2 ( i )-a22( i )«scal e 
20 continue 

cl ose (un i t=8 , status = ’ keep ' ) 



del=time(2)-timed) 
tt*de l*(n-l ) 

call fft(al .b,xl .nu.del ) 
call f ft (a2 .b, x2 .nu.del ) 

CALCULATION OF COHERENCE 

do 100 i=l.n/2 
rl =real (xl ( i ) ) 

6l-imag(xl ( i ) ) 
r2=real(x2( i )) 
s2= imag(x2 ( i ) ) 
pi 1 ( i )=atan( sl/rl ) 

p22(i)=atan(s2/r2) 
if (rl . It.O.and.sl .gt.O) 
if (rl . It.O.and.sl . It. 0) 
if (rl .gt . 0 . and . si . 1 t . 0 ) 
if (r2.1t.0.and.s2.gt.0) 
if (r2.1t.0.and.s2. It.O) 
if (r2 . gt . 0 . and. s2 . 1 t . 0) 


FUNCTION AND PHASE ANGLE. 


pi 1 ( i ) =pl 1 ( i )+ pi 
pH (i )=pll( i )+ pi 
pH(i)=pH(i)+ 2*pi 
p22( i )=p22( i )+ pi 
p22( i )=p22(i )+ pi 
p22(i)=p22(i)+ 2*pi 


tl=rler2+sl*s2 

t2=sl*r2-rl*s2 
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pcrm«=8qrt < t l**2 + t2«*2 > 
power! ( i ) ^epri ( rl««2-fel««2 > 
power2( i >*Bqrt<r2#«2+82**2> 
pcr< i >Bpcrm/ (power! ( i >*power2( i )) 

!00 continue 

do 50 i*2,n/2 

25 if (p! ! ( i > . gt . p! ! ( i -! > > go to 50 
p! ! < i )=p! 1 ( i )+2*pi 
go to 25 
50 , . . cont jnue 

do 60 i=2.n/2 

55 if (p22( i > .gt.p22(i-!>> go to 60 

p22< i )«p22( i >+2*pi 
go to 55 
60 continue 

do 65 i=2,n/2 
p!2<i)=<p!!<i)-p22<i)> 

65 continue 

open<unit = !7. f i le=’ax’ ) 
write<!7,555)<b(i>,p!2(i),i=2,n/2> 

555 format<5x. f!2.5.8x.f!2.5.4x, f!2. 5) 

close <unit=!7> 
stop 
end 

subroutine fft<a,b.x.nu.del> 
dimens ion x<!!00},b(!!00}.a<!!00> 
complex x.u.w.t 
n=2**nu 

pi=5. 141592655 
tt=<n-l)*del 

do 5 j j = l .n 

b( j j > = < j j-1 )/tt 

5 continue 

do 6 i=l.n 

x( i )=cmplx<a< i ) .0.0) 

6 continue 

Do 20 1=1 . nu 
le=2**<nu+l-l > 
lel=le/2 
u=(1..0.) 

w=cmplx<cos<pi/f loat< lel)>.-sin<pi /f loat< lei ) ) ) 


do 20 }=1 , lei 
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do 10 icj.n.l* 
ip= i + 1 el 
t*x( i )+x{ i p> 
x(ip>®(x<i)-x<ip>>#u 
10 x(i)*t 

20 u®u*w 

nv2*n/2 
nml =n- 1 

j = l 

do >0 1 = 1 ,nml 

if ( i .ge. j >go to 2? 

t=x( j > 

x< j )=x( i > 

x( i)=t 

25 k=nv2 

55 if(k.ge.j>go to 30 

j= j-k 
k=k/2 
go to 35 
30 j=i+k 

return 

end 

subroutine trap<x,n.e) 
dimension xdlOO) 

fx=0-0 

do 10 i=2,n-l 
fx=fx+x( i ) 

10 continue 

e=<x<l >+x<n) + 2.0*fx>/<2.0*(n-l )> 


return 

end 
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SAMPLE SEA- PROGRAMME TO GENERATE REGULAR WAVES 


^include <default.h> 

^include <defaul t . ttf > 

^include <macro8.h> 

run "1.00 hz wave 0.04m" with 
OneFreq__0<l .00,0.04) 

sample "Force__Transducer 7 April 1994" on 0 
sample "Wavegauge PI 7 April 1994" on 4 
end 


c 


Regular wave parameters 


frequency = 1.0 Hz 
amplitude = 4.0 cm 
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SAMPLE SEA- PROGRAMME TO GENERATE RANDOM WAVES (JONSWAP SPECTRUM) 


^include <default.h> 

#include <defaul t . ttf > 

^include <macro6.h> 

Mave Ml » )on6Map(0 . 7$ . 0 . 0081 , ^ . 0 , 0 . 07 . 0 . 09) 
wave w2 « random(wl , 1 > 

run "Jonswap spectrum" with 
wavemaker w2 with TTF on 0 

sample "Force_Transeducer 06 sept 1994" on 0 
sample "Wavegauge <p2> 06 sept 1994" on 4 
end 


c Oonswap parameters peak frequency = 0.75 Hz 

alpha = 0.0081 
gamma 5.0 
sigma a = 0.07 

sigma b = 0.09 



APPENDIX B 


Theoretical calculations for regular waves 


Exper i went 1 
Given 

f = 1.00 Hz 


Uaveheight H = 9.468 cm 

Depth of water d = 70.00 cm 

Dia of cylinder D = 5.00 cm 

Calculate, 


d 

70 

= 0.071 


981x1 

H 

9,468 

= 0.0097 

gT^ 

981x1 


Referring to the shore protection mannual, vol. II {19771 

L. 


1,00 for 


= 0.0097 (Figure 7.44) 


gT‘ 


where L. = Wavelength in feet for period T and depth d 
A 

2 

L = Deepwater wavelength <L = 5. 12 T ) in feet 
o o 

= 5. 12 X 1 = 5.12 

= 156 . 05 cm 

L . = 156 . 05 cm, 

A 


D_ 

L, 


5 


156.05 


k0.032 < 0.05 


Hence, Morison equation is valid. 


Now, 


Reynold 

U 

max 


number Re 


riH ^o 



U D 
max 

V 
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V 


Re 


TT X 9,468 
1 

kinematic vi 

2 . 

10 cm /s 

29.74 X 5 
10 ~^ 

14872 


29.74 cm/« 

scosity of water at 20 


O 


c 


Keu 1 egan-Carpenter number 
U T 


max 

D 

29.74 x 1 
5 

= 6.00 


In order to calculate maximum in-line force and maximum in-line 
moment. Dean’s stream function theory is being used here for Re = 
14872. 

4 

C. = 1.2 if Re < 4x10 according to Thirrot, et al 11971), 
o 

5 

C = 2.0 since Re < 2,5 x 10 according to Thirrot, et al (1971). 
m 



2,0 5 

1.2 ^ 9.468 

= 0.86 

Using Dean's curves given in SPM, vol. 11 (1977). 


for 




0.0097 

0.071 


and W 


0.88 



82 


Ue get 

4> = 0.3527 

m 

ck = 0. 2564 

m 

F = maximum In-line force 

m 

2 

= ^ w C . H D (w = unit weight of water) 

m a 

= 0.3527 K (1000 X 9.81) x 1.2 x <9,468 x lo"^ ) x 0.05 

= 1.863 N 

M = maximum In-line moment 

m 

= a w D d 
m d 

* 0.2564 X (1000 X 9.81) x 1.2 x (5 x lO"^) x 
X (9.468 X lO”^) X 0.7 

= 0,947 Nm 

Calculation of Lift force 

KeuI egan-Carpenter number 

K' = 6.00 

Using the curve for C given by Sarpkaya (1976a) 

L 

Lift coefficient C = 3,3 for Re = 14872 and K = 6,00 

L 

4 a = maximum lift force 
Lrl 

= . Pd D ^ 

2 

(O 2974) 

= 3.3 X (100 X 0.7 X 0.05) x 

* 5. 198 N 

In summary, 


(a) 

Max imum 

in-line 

force = 

1.863 

N 

(b) 

Max imam 

in-line 

moment = 

0.947 

Nm 


(c) Maximum transverse (lift) force 


5, 198 N 
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APPENDIX C 

Sbine Important def initions 


Wavehfeight : 

It is difference in elevation between the topmost crest 
and the bottommost trough of a wave over one wavelength. For a 
long wavetrain, waveheight is determined by averaging over 
individual wavelengths. 

Root mean square Crms!> value : 

The rms value of a variable x(t) which is function of an 
independent variable t is mathematically defined as, 

T 

X = vf x^<t).dt 

rms T J 

o 

where T (assumed to be large) is total t, over which we have a 
continuous record of x. 


Spectral width parameter C€> ; 


It is the measure of randomness of a random wave signal 
defined in terms of spectral moments as follows 


where m 

o, 

spectra. 

spectrum 



TOg are zeroth, second and 

e takes values between O and 1. 

€ 0 (Raleigh distribution) 


fourth moments of the 
For a narrow band 
whereas for Gaussian 


distribution e -t 1.0. 
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Spectral density t 

The spectral density function gives the relative 

importance of a particular frequency component in a random wave 

2 

signal. The area under the spectrum is the variance (<7 ) of the 
signal . 

Time period ; 

It is the time after which a signal repeats itself. 

Phase difference ; 

The angular separation between two harmonic waves or 
harmonic components of a random wave. 



